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We describe a study of the underlying event in jet and minimum bias events using 
data from the CDF detector. The underlying event contribution to the jet energy 
has been calculated in jet events and compared to the results of two Monte Carlo 
programs: Herwig and Pythia. The analysis has been carried out at two different 
center of mass energies: y/s =1800 and 630 GeV. For most observables, good 
agreement is obtained with at least one of the Monte Carlo programs. Neither 
program describes all features of the data. 



1 A measure of the underlying event in jet events 

In hadron-hadron collisions, in addition to the hard interaction that produces 
the jets in the final state, there is also an underlying event, originating mostly 
from soft spectator interactions. In some of the events, the spectator inter- 
action may be hard enough to produce soft jets {double parton scattering) U. 
This contribution may be calculated perturbatively. In order to compare the 
results with NLO QCD calculations, the contribution due to the underlying 
event must be subtracted from the jet cone. CDF assumes an uncertainty of 
30% on the underlying event contribution to jets. This is the dominant source 
of systematic uncertainty for the CDF inclusive jet cross section at low Et B. 

In order to estimate the contribution of the underlying event in jet events, 
we consider two cones of radius 0.7 at the same rapidity and at ±90° in 
azimuth from the most energetic jet in the events. The cone size is the same 
as used to reconstruct jets in the inclusive jet cross section analysis at CDF. 
We require the leading jet to be in the rapidity region \ri\ < 0.7. The two cones 
are located in a semi-quiet region, far away from the two leading jets, but still 
in the central rapidity region. For each event we label the cone which has 
the maximum ^pt {max cone) and the cone with minimum ^pt {min cone) , 
where X]Pt refers to the sum of the transverse momentum of the tracks inside 
the cone. The third parton, in NLO perturbative calculation of jet production, 
may fall in at most one of these two regions. In this approximation, the 
difference between the max and the min cone provides information on this 
contribution, while the min cone gives an indication of the level of underlying 
event. 
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Figure 1. On the left: 2_]P' i^^^ide the max and min cone in jet events at a center of 
mass energy of 1800 GeV as a function of Et of the leading jet. On the right: the Swiss 
Cheese configuration is shown as a function of the Et of the leading jet 



We use CDF data from Runlb, at a center of mass energy of 1800 GeV, 
and data collected in December 1995 from Runic at ^/s = 630 GeV. In jet 
events at 1800 GeV, we consider four different data samples, each with Et of 
the leading jet larger than 40, 75,|J.OO and 130 GeV. We use two Monte Carlo 
programs, Pythia Q and Herwig El, to generate four samples with the same 
cuts as the ones applied on data. The results from the Monte Carlo are then 
passed through QFL, the fast CDF detector simulation. We correct the tracks 
for the reconstruction efficiency and require pt track > 0.4 GeV/c. Systematic 
uncertainties have been estimated to be of the order of 10% or less. 

Figure ^ (left) shows X)Pt inside the max and min cone as a function 
of Ef of the leading jet. We use Herwig with its default parameters for the 
soft physics description and tune Pythia to attempt to reproduce the data. 
The regularizatioH scale of the transverse momentum spectrum for multiple 
interactions (pto) a has been set to 2 GeV/c. Herwig+QFL, Pythia6.115+QFL 
and the data have a similar behavior for the max and min cone: the min cone 
stays flat while the max cone increases with the Ej of the leading jet. The 
simulation reproduces the level of transverse momentum in the min cone in 
CDF data. In the max cone, Pythia is about 300 MeV larger than Herwig, 
with the data being in the middle. Since the min cone is well reproduced, this 
difference is probably due to the different algorithms used by the two Monte 



Carlo programs to generate parton shower. 

An experimental approach to study the underlying event, complementary 
to the max and min cones in jet events, is the measurement of the sum 
of the transverse momentum in the central region {\r]\ < 1) excluding the 
transverse momentum of the particles inside a cone 0.7 from the center of 
the two or three most energetic jets in the event. This configuration is called 
Swiss Cheese and is shown in figure (right). There is a good agreement 
between data and simulation when the three most energetic jets in the event 
are subtracted. Here, the contribution from a third parton in the event should 
be subtracted. Again, when only two jets are subtracted, Pythia is larger than 
Her wig. 

At ^/s — 630 GeV, the simulation agrees with data. Here, in order to re- 
produce the data, the regularization scale of the transverse momentum spec- 
trum for multiple interactions (pto) in Pythia has been set to 1.4 GeV/c. 
A dependence of pto on the center of mass energy has been implemented in 
versions of Pythia after 6.12. However, the dependence in Pythia is slower 
than what we found in this analysis and the default values underestimate the 
number of charged particles produced by multiple interactions. 

2 Properties of soft minimum bias events 

In CDF, the minimum bias trigger requires coincident hits in scintillator coun- 
ters located at 5.8 m from the interaction point. It samples double diffractive 
and non diffractive interactions from beam-beam collisions. The diffractive 
contribution to overall cross section is small. 

Minimum bias events have been analyzed at both 1800 and 630 GeV. 
Figure (left) shows the transverse momentum distribution of the tracks in 
the central region (|?7| < 0.7) at y/s — 1800 GeV. Pythia reproduces the CDF 
data well, while Herwig does not generate tracks with Pf > 5 GeV, indicating 
a lack of semi-hard physics description in its model for the minimum bias 
events. At y/s = 630 GeV, neither Herwig or the tuned Pythia reproduce the 
data correctly. 

We split the minimum bias sample in two subsamples: a soft one, with 
events without clusters of Et > 1.1 GcV in the rapidity region \ri\ < 2.4, and 
a hard one, with events containing at least one of such clusters. In figure ^ 
(right), the multiplicity in the central rapidity region {\ri\ < 1), defined as 
the number of selected tracks in the event, is shown, for the soft samples at 
1800 and 630 GeV, in KNO variables. The distribution shows the validity of 
KNO u scaling in soft minimum bias events, while the hard sample slightly 
violates the KNO scaling. 
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Figure 2. On the left: pt of the tracks in the central rapidity region (|r) < 0.7) at -y^ = 1800 
GeV in minimum bias events. The number of entries in the simulation has been normahzed 
to the number of entries in the data. On the top part of the figure, only the region with 
Pi < 1 GeV is shown in linear scale. On the right: multiplicity in the central rapidity region 
(Irjl < 1) in KNO variables. On the bottom of the plot the ratio of the the two distribution 
is shown. The two continuous lines delimit the band of all systematic uncertainties 



3 Conclusions 

An improved tuning of Monte Carlo programs and understanding of the un- 
derlying event and minimum bias events at the Tevatron, is also desired in 
order to study the complex physics environment at the LHC, where about 25 
minimum bias events are superimposed to any event of interest. 
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